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Abstract 

We construct a non-equilibrium steady state and calculate the corresponding current for a 
(^ [ mesoscopic Fermi system in the partition-free setting. To this end we study a small sample 

^N) . coupled to a finite number of semi-infinite leads. Initially, the whole system of quasi-free 

fermions is in a grand canonical equilibrium state. At i = we turn on a potential bias on the 
prt'. leads and let the system evolve. We study how the charge current behaves in time and how 

it stabilizes itself around a steady state value, which is given by a Landauer-type formula. 



> 



^ ■ 1 Introduction 

i-G . At the present time one can essentially distinguish two different ways of constructing non equilib- 

♦r^' riuni steady states (NESS) for composed systems. 

(^ . The first method consists of preparing a partitioned initial state for the total system containing 

"t^ ' several sub-systems, each of which being in a different state of thermal equilibrium, and then put 

them into contact with each other at t = 0, and let the coupled total system evolve in time until 
it reaches a steady state. In the mathematical physics community this method goes back to D. 
Ruelle |30] . |31j . It was seriously promoted during the recent years through numerous papers, see 
e.g. [291 Ql HSl [SB [m [1] and references therein. One can allow the carriers to interact in the 
sample [15], and the theory still works. Note that even if one chooses to turn on the coupling 
between the reservoirs in a time dependent way (for example adiabatically) , the results remain the 
same [T5] . 

CO ' The second method deals with those situations in which the initial state is an equilibrium state 

ly-v , for the already coupled (i.e. partition-free) total system. The partition-free approach goes back 

(^ ' at least to M.Cini |9]. This means that the initial state is not "partitioned" into a direct sum of 

^^ , equilibrium sub-states associated to e.g. different leads. The system is taken out of equilibrium 

by switching on an electrical bias between subsystems (leads) like for example turning on a d.c. 
battery, which in a certain way can be seen as changing the electro-chemical potentials of the 
)^j ■ leads coupled via a small sample. In contrast to the first method, there are almost no rigorous 

^ ^N , mathematical results on the second method beyond the linear response theory, or at least we are 

C^ ' not aware of the existence of such results. 

Although these two methods seem very similar, especially if one suddenly switches on the 
parameter bias in the partition- free system at t = 0, their implementations are different. One of 
the aims of this paper is to illustrate this observation. 

The main result of the present paper is that now we are able to construct a NESS and to study 
charge currents in the partition-free setting and for the full response. Let us describe in words 
what we do. 

For simplicity, in this paper we only consider two semi-infinite leads which arc both coupled 
with the same small sample when i < 0. The full system is in a Gibbs equilibrium state at a 
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given temperature and chemical potential. At i = we turn on a time-dependent potential bias 
V{t) between the leads, modeling a transient regime of a d.c. battery. At time ti > the bias is 
stabilized and remains constant in time afterwards. The statistical density matrix p{t) is found 
as the solution of a quantum Liouvillc equation, with an initial condition given by our global 
Gibbs state at t = 0. The time-dependent charge current from one lead to the other is defined 
as the mean value of a current operator in the state p{t), sec (j2.8|) - (|2.12|) for details. A priori the 
current depends on time, on the way we switch on the bias, and the point where we make the 
measurement. 

In Theorem 12.41 we show the existence and compute the ergodic limit of this charge current. 
The limit depends neither on the way we switch on the bias, nor on the point where we measure 
the current. We also obtain an explicit Landauer-type formula for this limiting charge current 
value, involving the transmission coefficients between the leads. 

Establishing Landauer-Biittiker type formulas (see e.g. [6j[7l[2ll3) starting from first principles 
but in the partition free setting was the original motivation of a number of remarkable physical 
papers, see for example [T3], [12] j [1] and references therein. Probably the state of the art of 
this subject seen from a physical perspective is to be found in two papers by Stefanucci and 
collaborators j33l I34j in which the partition free approach is combined with the Grccn-Kcldysh 
theory and a number of very current interesting formulas are proposed. 

A first mathematically sound derivation of the Landauer-Biittiker formula in the partition free 
approach under the linear response approximation was obtained in |llj and further investigated in 
[H] . In [TU] we significantly improved the method of proof of [TT] , which also allowed us to extend 
the results to the continuous case. Another challenging open problem is to extend the formalism 
in order to accommodate more efficient numerical current computations in transient regimes (see 
[24l [25l [26l [27] and references therein), and locally interacting fermions. 

The structure of the rest of the paper is the following: 

• In Section [5] we introduce the model and define the transient charge current in (|2.12p . The 
main result is formulated in Theorem l2.4l 

• Section [3] starts with a list of well-known facts about the spectral and scattering theory of 
mesoscopic systems coupled to semi-infinite leads. The second part of the section is dedicated 
to the proof of our main theorem. At the end we give a list of open problems. 

2 Set up and main results 

We work with a discrete model in a one-particle Hilbert space Ti. Following the physical convention 
we define the scalar product to be linear with respect to the second variable. 

Our carriers are quasi- free fermions (electrons). A small sample S is modeled by F C Z^, chosen 
to be a finite subset of Z^. We couple S to two "one-dimensional" semi-infinite discrete leads 
a = 1, 2. The sites of a lead (building its standard basis) are indexed by the set Nq, := {0, 1, 2, ...}. 
Thus, \ja) denotes the basis element at the site with the number j of the lead a, see Fig.l. The 
total one-particle Hilbert space is a direct sum of the space modeling the sample F C Z^, and two 
spaces corresponding to the leads {Nq,}q=i_2 : 

H:=l\r)®f{Ni)®l\N2) . (2.1) 

We denote by {\m,n)}(^jyi,n)£Z^ and by {|:/Q)}jQeNQ the corresponding orthonormal bases of the 
spaces /^(Z^) and Z^(Ma), where a = 1,2. 

Now we describe our one-particle Hamiltonian. For the sample S we may choose any self- 
adjoint bounded operator H^ . For example, we can choose H^ to be the restriction to Z^(F), of 
a lattice Harper- type operator with Dirichlet boundary conditions on F, but the concrete model 
for H^ does not play any role in the proof of our results. Notice that F is chosen to be finite, but 
can be arbitrarily large. 
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Figure 1: A sample F connected to two semi- infinite leads Ni and N2. Here \ja) denotes a basis 
element at the site with the number j of the lead a 



On each lead a = 1,2 we define the identical one-dimensional discrete Laplacians acting on 
the functions from /^(Nq.) with Dirichlet boundary conditions on Nq: 

{H^^){n):=tL{^{n + l) + ^{n-\)}, n > 1; (if^*)(0) := iL*(l), 






(2.2) 



where i^ > is a hopping constant. In the following, we denote the Hamiltonian corresponding 
to these three disconnected subsystems by: 

Hn := H^ + H^. 

The coupling between the sample and leads is described by the tunneling Hamiltonian (sec Fig.l): 



H^ := tE{|0„)(5"| + |5")(0„|} =: ff^^ + H 



SL 



(2.3) 



a=l 



Here r > is the hopping parameter between leads and the sample. The interaction (|2.3p simulates 
a quantum point constriction, or a tunneling barrier. Here |0q) is the first site on the lead a, and 
|iS") is the corresponding contact site \nia,na) on the sample coupled to the lead a. 
Then the total one-particle Hamiltonian takes the form: 



H ■.= H'^ + ^H^ + H^ = H^ +H^ + H^^ + H^^=: Hq + H'^ 



(2.4) 



Remark 2.1. As we mentioned before, our results can be extended through verbatim to more 
general choices of Hamiltonians H^ and TJ"^. The key properties that we need are: the absolutely 
continuous spectrum in the leads, and a finite rank operator coupling between a finite sample and 
the leads. 

2.1 The state and charge current 

At t < the total coupled system (|2.4[) is at equilibrium for a given temperature 1//3 > and a 
chemical potential /x. Since we work with non-interacting fermions, the corresponding one-particle 
Fermi-Dirac equilibrium density matrix is the operator 

/W= e/3(H-l) + i '^^^' (2.5) 

defined on the Hilbert space (|2.ip . 

At the moment t = we turn on a bias on lead number one in the following way. We fix <i > 
and choose a real and continuous function (j) which has the property that 0(t) = if i < and 
(/>(t) = 1 if t > ii. Let u > 0. Denote by Pi : H i— > '^(Ni) the projection on the lead number one. 
Then define the time dependent potential bias as: 

Vi{t) := v^Pi . (2.6) 

Denote by U{t) the unitary evolution associated to H+Vi{t) through the time dependent Schrodinger 
equation: 

idtU{t) ^{H + Vi{t))U{t), t/(0) = I . (2.7) 

The density matrix at time t > is a solution of the Liouville equation and be expressed by: 

Pit) :- U{t)f{H)U{tr . (2.8) 

Denote by Pj " ■ ^ ^^ ^'^{{n, ri + 1, . . .}), the projection on the second lead from which we exclude 
the first n sites. If n = 0, then it is just the projection P2 on the lead 2. We define the current 
operator modeling the measurement of the charge flow at site n by: 

j„ :=z[i/ + yi(i),Pi"'] =i[i7,Pi"^] , JO = *[i?^,P2] . (2.9) 

Remark 2.2. Clearly, the current operator has finite rank, thus it is trace class. This is one 
important feature which is only true in the discrete setting. It significantly simplifies the technical 
estimates compared to the continuous case. 

Remark 2.3. To obtain (|2.9p we used some evident support properties of the projections P2 , 
which imply: 

[Pi,P2^")]=0, P2("^P2("+^U P2^"+^\ [iJ^,P2<"'] = 0, Vn>0, (2.10) 

and also: 

[H^,P2\ = , [H^, P^"^] ^ 0, [H'^, Pi"^] = V?7, > 1. (2.11) 

These properties make the commutators in (|2.9p nontrivial, and are important for the study of 
the current propagation (see the proof of the point (ii) of our main theorem) . 

The charge current flowing through the second lead at time i > and measured at site n is 
the expectation of the operator j„ from (|2.9p in the quasi-free state defined by the time dependent 
one-particle density matrix p(t): 

/(<,n):=Tr{p(t)j„}. (2.12) 



2.2 The main theorem 

Now we are ready to formulate our main results, collected in one theorem. 

Theorem 2.4. (i) The following ergodic limit exists and is independent of n, ti and <j): 



1 '■^ 



loo ■■= lim - / I{t,n)dt. (2.13) 

(ii) Fix t > 0. Then the current vanishes if we measure it infinitely far inside the lead 2: 



lim /(t,n) =0. (2.14) 

(iii) Assume that the operator H + vPi has only finitely many eigenvalues. If we measure the 
current very far inside the second lead (hut not infinitely far), the transient/ oscillatory effects will 
he weaker and weaker and the current defined in 112. 12\) will slightly fluctuate around the value 
112.18]) . More precisely: 

lim limsup |/(t, n) — /oo| = 0. (2-15) 

n->oo t-^oo 

(iv) Denote hy T^v" W ^^^ transmittance coefficient hetween the two leads at hias v for the spectral 
parameter X [see (j3.19p for a rigorous definition). Then we estahlish the following Landauer-type 
formula: 

Ioo=2tt [ {/(A) - /(A - v)}rl^\x)dX . (2.16) 

Remark 2.5. From (ii) one concludes that if we fix the time of measurement t and push the 
measuring point n to infinity, the current tends to zero. Although this does not imply that the 
propagation speed of the current is finite, a property which we cannot expect to hold true because 
of the non-relativistic dynamics on the leads. On the other hand, in (iii) we prove that if the 
current measuring device is placed further and further away from the sample, after waiting a 
very long time the current becomes non-zero and has weaker and weaker fluctuations around its 
steady-state mean value. In other words, the limits t — ;> oo and n — ;■ oo do not commute. 

Conjecture 2.6. As a complement to Remark \ 2. 5[ we conjecture that the group-velocity of the 
spatial correlations in our model is finite, i.e. if A is any ohservahle supported in a neighhorhood 
of the sample, then there exist some positive constants C , M and v such that the Lieh-Rohinson 
type hound j23j : 

II [A, U{typ!^''^U{t)] II < Ce-^^f"-"*' (2.17) 

holds true for every t > and n > 0. The exponential hound sounds as a strong one, i.e. one can 
not exclude a priori a power-like decay. 

Remark 2.7. Our proofs are exclusively based on one-body scattering methods. We do not use 
the many-body language, which is unavoidable only if the carriers interact. 

Remark 2.8. By the same reasons, one can completely characterize the many-body states aj(-) 
on the Fermi algebra CAR(7^), (algebra of the Canonical Anticommutation Relations) by a one- 
particle density- matrix operator p defined otiT-L. If at i = the state w(-) is the grand-canonical 
equilibrium state on CAR(H) of a non-interacting Fermi system (|2.4I) (equilibrium quasi-free 
state), then the density-matrix operator is equal to p{t = 0) = f{H) (|2.5p . The evolution (j2.7p 
preserves this property, i.e. it transforms this state into a non-equilibrium quasi- free state (J2.8I) . 

3 Proof of the main theorem 

We start this section with a list of well known facts about the spectral and scattering theory of 
mesoscopic, systems coupled to semi-infinite leads. This will help us to fix notation and streamline 
the proof of the theorem. 



3.1 Some spectral and scattering background 

First we recall some elements of the stationary scattering problem associated with the pair of 
Hamiltonians {H + vPi,Ho + vPi), where H := Hq + H^ . In this case the free system consists of 
the leads with a bias v localized on the first lead together with the decoupled inner sample, and it 
is described by the Hamiltonian H^ + vPi. The perturbed system also contains the coupling H^ . 

The operator i/o = iJ^ + if'^ has as a subspace of absolute continuity H''^(ifo) = ®a=i^^(Na)- 
Since the operator H — Hq is of finite rank, the trace class scattering theory implies that the M0ller 
wave operators 

f]^"^ = 5-lim e''*(^+^^i)e-**(^«+^^i)^ac(i?o + vP^) , (3.1) 

exist and are complete, see e.g. [32], [3n|- Here E^c{Ho + vPi) = -E'ac(-ffo) denotes the projection 
on the absolutely continuous subspace 'H°''^{Hq + vPi) C %, or 'H°''^{Ho) CH oi the corresponding 
operators. The location and nature of the spectrum of operators like H was extensively studied in 
[llj : one can prove under generic conditions that there arc only finitely many eigenvalues, while 
the singular continuous spectrum is always absent. 

It is known that the set of (normalized) generalized eigenfunctions of H^ on the semi-infinite 
leads a = 1,2 have the form: 

vI'.(A)=^*(A;m)|m„), *(A; m) = ^^^^S^^^. (3.2) 

^0 V^tL sm(fc) 

Here the spectral parameter A = Xk{'-= 2iiC0s(fc)) for k G (— 7r,7r). The generalized Fourier 
transformation associated to these eigenvectors is defined by 

2 2 

F: ^l'{N^)^^L'-{[-2tL,2tL]), (3.3) 

a— 1 a— 1 

[F(<i>)]«(A) = (vI/,(A), $«)P(N„) - Y. WM^aim). (3.4) 

Its adjoint is given by 

F* : L\[-2tL, 2tL]) ^ Pin^), (3.5) 

a— 1 a—1 

f2tL 

[F*{E)Um)= S„(A)^(A;m)dA. (3.6) 

J-2tL 

We see that F is a unitary operator, and that FH^F* is just the multiplication by AI on the space 
which is a direct integral L_2f 2t 1 ^^'^^ ~ ®q=i ^^([~2tL, 2<l]), i-e. 



m>0 



2 2 



FH^F* = XI dX. (3.7) 

J[~2tL.2tL] 

If the bias is present on the first lead, the situation is changed. Since H°''^{Ho) = 'H'"^(i?^), the 
generalized eigenfunctions of Hq + vPi are chosen to be 

*J''^(A;?«) :=*(A-w;m), X e [-2tL + v,2tL + v], m > 0, (3.8) 

^^''^(A;m) :=*(A;m), X e [-2tL,2tL], m > 0. (3.9) 

The corresponding generalized Fourier transformations are: 

2 

F, : l^Na) -^ L^{[-2tL + V, 2tL + v]) © L\[^2tL,2tL]) (3.10) 

[F,($)]i(A) - J2 *(A-v;m)$i(m), [F„($)]2(A) = ^ *(A;to)$2(™) ■ (3.11) 

m>0 m>0 



Therefore, we can construct generalized eigenfunctions oi H + vPi , as solutions of the Lippniann- 
Schwinger equation: 

-^i^H^; •) - ^L^'HA; •) - (ffo +vP,-X- iO+)-'H^-^i-\X; •) . (3.12) 

These generalized eigenfunctions have the following very useful intertwining properties between 
the subspaces of absolute continuity of the operators Hq + vPi and H + vPi , which can be formally 
written as: 

$(:nA;-)=^V"'*H(A;.), (3.13) 

¥:\Xr) = {n^_^^*^^:\X;-). (3.14) 

The scattering operator S*^"^ : 'H°-^{Hq + vPi) i~> H'^'^{Hq + vPi) is a unitary map acting 
on n'^'^iHo + vPi) = W'^iHo) = 0^^i P(N„), and it is given by S^'''^ = {f]^''^}*^!^''^. Then the 
corresponding transition T-operator is defined by 2'KiT^'"^ := I— S**-"-* . In the spectral representation 
of H^ + vPi in the space Ji_2t 2t 1 ^^^-^j the T-operator is a A-dependent 2x2 matrix with 

elements denoted by ^^^(A). Using p.6p one gets the representation: 

Y^ iS(A)S^(A) = -i.[F(I- 5(^'))F*S]„(A) . (3.15) 

13=1,2 

Then with the help of the generalized eigenfunctions we can express the T-matrix elements as: 

iS(A):=(*i''^(A;-),i?'^<i>^"'(A;-))- (3.16) 

Since S is unitary, one gets the relation i{T -T*) = 2ttT*T = 2'kTT* (Optical Theorem), which 
implies: 

Im{4^,)(A)} = ^ (|42^(A)P + |42Ha)P) , (3.17) 

l4i^(A)P = l42'(A)p. (3.18) 

The transmittance T^J (A) between the leads a and /3 for a given energy A is defined by: 

'rj;^(A):=|iS(A)r (3.19) 

Note that by definitions ^^, dSSj) and ((3?T6| the transmittance 7^2''' (A) = if A ^ [-2*^ + 

v,2tL+v]r\[-2tL,2tL\. 

3.2 Proof of (i) 

By (|2.6I) and (|2.7p the evolution operator U{t) obeys for t > ti the equation: 

U{t) == e-''(*-*i»-^+"-Pi'[/(ti) . (3.20) 

Then by H == W"p{H + vPi)®W''{H + vPi) and by (E^ we obtain for the current ((^1^ measured 
at site n the representation: 

/(i,n) = Tr{e-*(*-*i)(^+^-f'i)j7(ti)/(iJ)t/*(ii)e*(*-*i)(-f^+"-f'i\7„} 

= Tr{e-*(*-*i)(^+"^i)t/(ii)/(i7)f7*(ti)e*(*-*i)(^+"^i)-Bpp(i7 + wPi)j„} 

+ Tr{e-*(*-*i»^+"-f'i)[/(ii)/(i7)t/*(ti)e^(*-*i)(^+^-f'i)-Bac(ii" + vPi)j^} 
^■.Ipp{t,n)+h,{t,n), (3.21) 



where Epp{H + vPi) denotes the projection on the pure point subspace WP{H + vPi). By virtue 
of (|2.9p one gets the identity: 



e 



,(t_t,)(ff+„P,)_g^^^^ + «Pi)jne-^(*-*i)(^+"^i^ 



PPV 

'^ {£;pp(i7 + vPi)e^(*-*i)(^+"^i)P2^"'e-^(*-*i'(^+^^i)}. (3.22) 



dt 

Let us for now assume that H + vPi has a finite number of eigenvalues. This means that Epp{H - 
vPi) is trace class. Now if T > ti, the pure point part of p.2ip yields: 

T rti rT rti 



Ipp{t,n)dt= I Ipp{t,n)dt+ / Ipp{t,n)dt = I Ipp{t,n)dt 
Jq Jti Jo 

+ i:Y{U{t{)!{H)U*{U)D{T)Epp{H + vPi)} , (3.23) 



where the operator D{T) := e'(^~*i)(^+^'f'i)P2^"^e-*(^-*i)(^+"^i) - Pj^"^ is uniformly bounded in 
T. Since the first integral in the right-hand side of p.23[) is finite, the pure point spectrum does 
not contribute to the ergodic limit p.l3p . In the case when Epp{H + vPi) does not have finite 
rank, we have to employ an e/2 argument based on the fact that Epp{H + vPi)jn can be arbitrarily 
well approximated in the trace norm with an operator containing the projection on a sufficiently 
large (but finite) number of eigenvalues of _ff + vPi . This approximation will be independent of 
T, so the previous argument can be repeated. 

So, it remains to investigate /ac(t, n) and to show that it actually converges when t -^ cxd. To 
this end we start with three technical lemmas: 

Lemma 3.1. The operators U{ti) - e'^HoU-tvP, /Ji 0(r)dr ^^^ U*{ti) - 0""''^^+''"^^ ^«' "^(^)''^ are 
compact. 

Proof. Since the following Dyson-type equation: 
is equivalent to 

Uiti) = Q-^Hoti-ivPi /qI ^{T)dT _ I I ^-iHo(ti-t)-ivPi /^ i 't>{'r)dT jjT jj f.\^^ 



Jo 

we use that i?^ is a compact (finite-rank) operator in order to finish the proof. D 

Lemma 3.2. The operator U{ti)f{H)U*{ti) — f{Ho) is compact. 

Proof. It is an easy consequence of Lemma [XTl of the fact that Hq commutes with Pi, and of the 
observation that the difference f{H) — /(i?o) is a compact (even trace-class) operator. D 

Lemma 3.3. Let K be a compact operator. Then the following trace-norm tends to zero: 

hm ||i^e*(*-*i)(^+^'^i)i?ac(ff + t'Pi).?„||i = 0. (3.24) 

Proof. Since jn is from the trace-class (finite rank in our case), by standard e/2 arguments we 
can assume that operator K has a finite rank. Then the proof is a consequence of the Riemann- 
Lebesgue lemma. D 

Corollary 3.4. Use the identity U{ti)f{H)U*{ti) = {U{ti)f{H)U*{ti) - /(i?o)) + f{Ho) in the 
representation of Ig^{t,n). Then Lemma \3.2\ and Lemma \3.S\ imvlv the limit: 

hm \h,it,n) - Tr{e-'(*-*i)(^+"^i)/(iJo)e*^*-*^)(^+"^^)i?ac(i? + vPi)]n}\ = 0. (3.25) 



Now, to prove the ergodic limit (j2.13p it is enough to check that the trace appearing in p.25[) 
converges when i — >■ cx). To this end we use a standard trick of inserting the free evolution and 
then to use the identity: 

in p.25p . Using (|3.ip together with the fact that the wave operators are complete thus unitary, 
we obtain the existence of the following strong limit: 

where the limit operator projects onto H^^XHo). Finally, because jn is trace class we can conclude 
that the limit: 

lim 4c(t, n) = Tr |o^'V(-ftf^){^i''^}*-Bac(-ff + vPi)jn} , (3.26) 

exists and is finite. 

Remark 3.5. We were able to replace f{Ho) = f{H^) © f{H^) by f{H^) because the inner 
sample is projected out by the wave operator {n_^ }* on the right. 

Until now we proved that the ergodic limit p.l3p is independent of 4> and ti. The independence 
of n follows from the next lemma: 

Lemma 3.6. For any n > 1 one can establish the following continuity equation: 

Tr {n^^^f{H'^){n'-_;:'>}*E,,{H + vP,)jr,} = Tr {fi'^'V(i7^){^'^')}*i?ac(i? + vP^)Jo} 

-: /oo. (3.27) 

(n) 

Proof. Denote by Xn '■— P2 — P2 the projection on the first n sites of the second lead. Then 



(j2.9p and Remark 12.31 vield jo — jn = i[H + vPi,Xn]- Hence, the identity p.27p is equivalent to 
Tr {f]V"'/(i?'^){^+^}*^ac(i/ + vPi)[H + vPi,Xn]} = 0. 

But the operator Xn is trace-class, so we can undo the commutator. The wave operators intertwine 
between H+vPi and Hq+vPi, and Hq+vPi commutes with H^. It follows that H+vPi commutes 
with r2^''/(-ff^){^+^}*£'ac(-ff + wA)- Then the trace cyclicity finishes the proof of the lemma. D 

3.3 Proof of (ii) 

We start by proving that the current at t = (i.e. at equilibrium) is zero for all n. Indeed, 
according to (EH) and (1^1^ one has /(0,n) = Tr(/(7?)j„) = iTT{f{H)[H, P^""^]) for ah n > 0. 
Now fix n and denote by Hn the Dirichlet restriction of the operator H to the finite leads of 
length A^ < 00, where n < N. Denote by Pj the finite-rank projection corresponding to the 

(n) 

restriction of ij to the bounded second lead. Then one can prove a certain "thermodynamic 
limit" result [II]: 

/(0,7i)= lim Tr\f{HN)i[HN,Pt''''']}- (3.28) 

To understand why p.28p holds true, note that j„ = i[HN,P2 ] = i[-^>-P2 ] i^ ^ finite-rank 
operator which is independent of A^. Moreover, f{H) and f{HN) differ significantly from each 
other only very far from the support of j„- Details can be found in |llj . 

But Tr lf{HN)i[HN, ^2^"^'^] | = for aU A^ by trace cychcity. Thus fX^ shows that /(O, n) = 
for any n > 0. 

The next step of the proof is to show that for alH > one has: 

lim |/(t,?i)-/(0,n)| = . (3.29) 



Then by /(O, n) = for all n, the limit in p.29p would imply (j2.14[) . First we present the difference 
in ((3:291 as: 

Iit,n) - 7(0, n) = Tr{C/(t)/(ff )[[/*(<) - e'*^]j4 + Tr{[C/(f) - e'^'"]fiH)e^'"j„}. (3.30) 

Then we express the propagator U{t) (12.71) with the help of its time-ordered Dyson series: 



U{t) = e 



-itH 



(3.31) 



E 

fc>i 



'0 Jo 



k\ 



drfeT{0(Ti)e'"^^Pie 



—iriH 



.4{Tk)e'^''" Pie 



^-ITkH^ 



whereas its adjoint is given by: 



(3.32) 



E' 

fc>i 



fc! 



dri / dT2. 
Jo 



dTfeT{<?!)(Ti)e*^i^Fie''^i^...0(Tfc)e'^'=-f^Pie-^^^-f^}e"-f^ 



where T means time-ordering in decreasing order and T means time-ordering in increasing order. 
Note that in the formula p.30p the first term on the right-hand side contains the operator 



[U*{t) 



itH 



]jn', we want to show that its trace-norm goes to zero with n. By a simple support 



property (|2.10p one has P 

itH] p("-l) 



(n-1) 



jn = jn and since j„ is a finite-rank, it is enough to prove that 
[C/*(i) — e^'^jFj" ''' converges to zero with n in the operator norm. To this end we need a technical 
estimate given by the following lemma: 

Lemma 3.7. For any fixed t > one has: 



lim sup 

"^°°|r|<t 



irH p(n-l) 



Pie"" P. 



Proof. Since the operator H is bounded, for any e > there exists N^ such that 



sup 

\r\<t 



.rH_J2 



' i^T^H^ 



k=0 



k\ 



< e. 



(3.33) 



(3.34) 



The support properties (Remark [23]) and the one-step hopping in the Hamiltonian H imply that 
PiH''P^"^^^ =Oiin> N,>k. Hence, by (|53i)) we obtain that for n > N^ 



sup 

kl<t 



irH p("-l) 



Pie""P 



<e, 



which proves the lemma. 



(3.35) 
D 



itH] p("-l) 



Applying this result to the expansion (|3.32l) . one finds that [U*{t) — e^'^jFj" "' converges to 
zero in norm. This convergence allows to bound from above the limit of the difference (j3.30p 



limsup|/(t,n) -/(0,n)| < limsup |Tr{[t/(t) 



-itH 



]fiH)e^'"jn} 



(3.36) 



To estimate the limit p.36p we use the representation [U{t) 



-itHl 



-itH] 



{{^-P, 



(["/2])n 



/(i/)e**«j„ = [U{t) 



P2 }/(-ff)6'*^Jn- Since again the function /(_ff)e'*^ can be approximated 



in operator norm by polynomials in H, we can apply to (I — P^ )f{H)e.'^^^jn the same line of 

reasoning as in Lemma [3.71 to establish: 



lim 11(1- P2(["/'lV(^)e'*''i^i""'^ 



0, 



(3.37) 



([n/2]) 



j(»-l) 



since the distance between the supports of I — Pj and of P2 tends to infinity. For the term 

[U{t) - e-'*"]P^^"^^^'^ f{H)e'*" jn we use the representation ([33T|) and Lemma [3Jl which imply 



that the norm of [U{t) — e 



~itH]piln/2]) 



goes to zero with n. Together with p.37p this proves that 



the limit of the right hand side of p.36p equals zero, thus p.29p follows. 
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3.4 Proof of (iii) 

By virtue of p.26p and p.27p one has Mvat^oo Iac{t,n) = loo- Therefore, it only remains to 
estimate the current Ipp{t,n). This gives by p.2ip : 

sup|/pp(t,n)| < \\Epp{H + vPi)jn\\i . 

t>0 

Note the right-hand side of this estimate can be made arbitrarily small by increasing n, since we 
assumed that we have finitely many eigenfunctions which are necessarily localized near the sample 
S, thus 

lim \\Epp{H + vPi)P^''~^^\\^0. 

n— ^oo 

This finishes the proof of (iii). Note that in the exceptional case in which H + vPi could have 
infinitely many eigenvalues, this argument fails. 

3.5 Proof of (iv) 

To calculate the steady charge current (|2.16p we use our main formula (j3.27p in the form: 

I^ = Tr{r!V"V(^'^){"V"'}*^ac(i/ + ^^Pi)jo} = Tr{/(if^){r!!;'}*ii;ac(i/ + ^;Pi)jor!^"^} • 

Now using the spectral representation for if^ + vPi one can evaluate the trace on Z^(Ni) © 
Z^(N2) with the help of its generalized eigenfunctions (see p.8p and p.9p ). Then we obtain the 
representation: 

2tL+V , 

dX /(A - v) {^[^\X; •), {n^:^*E,,{H + vP,)jon^^^^[^\X; •) 
dX /(A) Ui^\X; •), {d^^*E,,{H + vP,)jon^:^^i^\X; •)) • (3.38) 

By (j3.13p for the scalar product in the first integral we get: 

^(¥,^\X■,■),{^^:^r[H'^,P2]n^:^^['\X■,■))=2Im(<^[^\x■,.),P,H^<^[^\x■, 
= 2lm(^P2{^[^\X;-)-{Ho+vPi-X-iO+)-'H^<^^^\X;-)},H^<i>^^'\X;-)) , (3.39) 

where in the second line we used the Lippmann-Schwingcr equation p.l2p . Notice that the vector 
P2{Ho + vPi - X- iO+)-^H^'^[''\X;-) S ^(Na) for almost every A. By Remark O one has: 
[Ho, Pa] = and P2P1 = 0, which implies P2 {Ha + vPi-X- iO+)-^H'^ =P2{Ho-X- iO+)-^H^. 
Taking this and the identity: Pa^j (A; ■) = into account, we can use the spectral representa- 
tion of Hq and decomposition the vector P2{Hq — A — iO+)~^iJ"^$J^" (A; •) over the generalised 

-2 (•^'; ■)}A'e[-2ti,,2ti,] 



eigenvectors {*^"'(A'; ■)}A'ef-2tz,,2tz,l to obtain 



2Im(P2{*(")(A;-) - {H^ + vP, - A - iO+)-iH^$l'')(A; •)}, P^^$("'(A; ■ 

= -2,rT^,->(A)x|-2i,,2.,lW- (3-«) 

For the last equality we used the Sokhotskii-Plemclj formula, and definitions (|3.16p . p.l9p . 
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By the same line of reasoning one gets for the second integrand in p.38[) : 

= 2Im{42'(A)} - 2Im^P2(i^o - A - iO+)-^H^'^i''\X; ■),ff^4"'(A; 

= 21m{il\x)} - 2n\il\x)\' = 2nrl^\X) , (3.41) 

where for the last identity we used p.l6p and p.l7|) . Note that here A e [— 2ti, 2ii]. Taking into 
account the symmetry ((XTOl) and plugging ([3^, ((XiT|) into (jX^ . we obtain ^JE^. 

Recall that r/a^^lA) = 0, if A ^ [~2tL,2tL] D [-2tL + v, 2tL + v]. D 

4 Concluding remarks 

In the present paper, we established a Landauer-type formula for the stationary current running 
through a discrete system with a (small) sample coupled to one-dimensional infinite leads. We 
give a rigorous proof of the existence of the ergodic limit of the charge current and then its explicit 
expression. Our strategy is based on the partition-free approach, it is quite general and demands 
a minimal information about the sample. 

There are several open problems which deserve to be mentioned. 

1. One of them is our Conjecture 12.61 about the Lieb- Robinson type correlation group velocity 
bound, which up to our knowledge it has not been studied before in this context. 

2. If V{t) is a time dependent bias between which after t = ti becomes a perfect monochromatic 
signal like Vq-I-Vi cos(aj<), then the ergodic limit exists and is independent of ii and of the site 
where one measures the current. Furthermore, the ergodic limit is given by a Landauer-like 
formula [18j . 

3. A computation of the current l{t,0) (see (|2.12p ). by expressing the evolution unitaries 
through the functional calculus associated to the resolvents, and the resolvents with the 
help of the Feshbach formula as in [TTl [T2]. Can one obtain an "easy" formula for I{t, n) at 
a given n? Can one study numerically the transient effects and check point (ii) in Theorem 

mp 

4. Study the resonant transport in the case of small coupling (0 < t << 1 and v a variable 
parameter). 

5. What happens with point (ii) of our theorem if there are infinitely many eigenvalues? 

6. Study the "wide band limit", or t^ — s- oo. 

7. Compute the first few corrections in v of the conductivity tensor. 

8. Introduce a Kohn-Sham interaction in the sample, as in Stefanucci's papers [12]. How can 
one properly formulate the mathematical problem in this non-linear case? Can one still 
prove the existence of a steady state? Is it unique? 
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